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I. INTRODUCTION 



Singlet states of the helium atom, especially its ground state, are best suited for precision 
studies of the electron-electron interaction at low energies. Recent measurements of l 1 ^ — 
2 1 P UJ and l 1 ^ — 2 X S [f2] intervals in helium reached the precision of about 10 ppb. Helium 
ground state ionization potential (the difference between ground state energies of the singly 
charged ion and of the atom) extracted from those measurements constitutes 



15-2P, 



^exp 

and 



l 1 ^) = 5 945 204 238 (45) MHz (1) 



15-25/ 
exp 



-(l 1 ^) = 5 945 204 356 (48) MHz, (2) 



respectively. 

Theoretically, the ionization potential can be calculated as a power series in the fine 
structure constant a. Leading 0(ma 2 ) contribution to the ground state energy is the lowest 
eigenvalue E of the nonrelativistic Hamiltonian 

entering into the Schrodinger equation Hip = Eip. Here and below I use the following 
notations: r a and p a are the position and momentum operators for the particle a with 
mass m a and electric charge z a (in units of the proton charge). The relative position of two 
particles is r a i, = r a — r^; for any vector v, v denotes \v\. The helium atom consists of two 
electrons with masses m\ = m 2 = m and charges z\ — z 2 — — 1, and the nucleus with mass 
m 3 = M and charge z 3 = 2. A general case of z 3 = Z takes into consideration helium-like 
ions with Z < 1/a. In the center-of-mass frame, total momentum of the atom vanishes, 
Y^aPa = so that only two of three position vectors are independent: I] a m a r a = 0. In 
singlet states, spins of the electrons sum up to zero while the orbital part of the wave function 
is symmetric with respect to permutation of the electrons positions, ip{ri, r 2 ) = ^(r 2 , ri). 

Relativistic and radiative effects shift the Schrodinger value of the energy by corrections 
of higher orders in a. In particular, the leading correction, 

SWz^/^jl z J(r 13 ) + 8(r 23 ) , na6(r 12 ) 



m\ 2m 2 m 2 

a ^ z a z h 



2 a - b>a m a m b 



r ab r 6 ab 



(4) 



arises as the average value of the Breit perturbation (see, e.g., 0) over the nonrelativistic 
wave function and is of the order a 2 relative to the Schrodinger energy. It is taken into 
account in (^) that the total spin of electrons and the spin of the nucleus are both equal to 
zero. 

A particular class of corrections appears due to the nucleus structure. There, the most 
important (and sufficient to be included at the present level of accuracy) is the effect of the 
nucleus charge radius i?N, 
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w? = ^i4E(<^N)>. (5) 

6 a 

Here r a N denotes the position of the a th electron with respect to the nucleus. 
The most recent theoretical result for the helium ionization potential, 

u^S) = 5 945 204 226 (91) MHz, (6) 

obtained in |4| includes along with @ and (|5|) the ma 5 order, the leading part of the ma 6 
order, and some estimates of higher order contributions. Its uncertainty is twice as large as 
that of the experimental data ([!]) and (0). The main source of the uncertainty in is the 
yet uncalculated part of the 0(ma 6 ) correction. 

The present work is the first of two devoted to the calculation of the helium ionization 
potential with 0(ma 6 ) accuracy. It contains the analytic part of the calculation and its 
main result is a set of effective operators which produce all 0(ma 6 ) corrections to singlet S 
levels of helium atom and low-Z helium-like ions. To make the presentation self-contained, 
I also briefly outline how the lower order correction can be obtained by the same method. 
The second paper || contains numerical results for the average values of the derived effective 
operators as well as all other known contributions to the helium ionization potential. 

The rest of this paper is organized as follows. Section describes general features of 
the approach. Order ma 5 effective operators are derived in Sec. |HI[ Sections |V| and [V] are 



devoted to the 0(ma 6 ) effective operators appearing from hard and soft scales, respectively. 
The final result of this paper is presented in the Conclusion. 

II. FRAMEWORK OF THE CALCULATION 

Since the early days of quantum electrodynamics (QED) the nonrelativistic expansion 
of an atom's ground state energy is known to break down at the ma 5 order ||. In contrast 
to the 0{ma A ) effective operators whose average values (|) are completely determined by 
the soft (p ~ 1/r ~ ma) scale, the operators of the next order in momentum-to-mass ratio 
are too singular to ensure finiteness of their average values over the ground state described 
by the wave function iftf\. It means that those operators become sensitive also to the hard 
(p ~ m) scale which is beyond the scope of the nonrelativistic expansion. Another novel 
feature of the O (ma 5 ) contribution to the energy is that the very picture of interaction 
between particles through a potential fails: virtual transitions from the atom's ground state 
to excited states and a photon become relevant. Thus, one more scale comes into play - this 
intermediate ultrasoft photon has an order ma 2 energy. The most natural way to calculate 
such a multi-scale shift of the energy is to divide it into several pieces each originating from 
its own scale and then use simplifying approximations suitable to that scale. For example, 
the nonrelativistic expansion is applicable at the soft scale. Alternatively, one can neglect 
bound-state effects at the hard scale. If all relevant contributions are included, their sum 



See Section III for details. 
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is independent of the details of the division. Applied to bound state problems in QED, 
this idea was first formulated in Ref. |7j as the nonrelativistic quantum electrodynamics 
(NRQED). 

Traditionally, in atomic calculations involving several scales some auxiliary parameter 
is introduced to separate a contribution of the given scale from the others (see, e.g., Ref. 
H where such a scheme is applied to the helium problem). For example, to divide soft 
and hard scale contributions one introduces A satisfying ma cACm and gets a final A- 
independent energy shift as a result of cancellation between two A-dependent contributions. 
The soft scale one includes A as the ultraviolet cutoff which makes average values of singular 
operators finite. Simultaneously, A cuts off otherwise infrared divergent on-shell scattering 
amplitudes which represent a hard scale contribution. The higher is the order of a, the more 
severe are singularities of both contributions and the larger is a number of A-dependent terms 
to be canceled in order to get a final result. The problem seems even less tractable when 
the wave function is known only numerically. 

Precise calculations of the positronium spectrum PHTSf have shown that contributions 
from various scales can be separated much more effectively by shifting the number of spatial 
dimensions d from three, d — > 3 — 2e. For consistency, the number of space-time dimensions 
in hard scale calculations should be shifted from four to 4 — 2e. This shift implies essentially 
that all objects defined originally for d = 3 are analytically continued to the complex plane 
of d. The main advantage of the dimensional regularization over the traditional scheme 
is that due to the analytic continuation all power divergences automatically drop out of 
calculations and one only has to keep track of logarithmic divergences which show up as 
inverse powers of e. 

Recall that there are two kinds of effective operators in NRQED. Operators coming from 
the hard scale are contact, i.e., they are proportional to delta-functions of distances between 
particles. Infrared divergences typical in hard scale contributions manifest themselves as 
inverse powers of e in coefficients of those delta-functions. On the other hand, the soft scale 
effective operators have finite coefficients at d — > 3. Ultraviolet divergences inherent to soft 
scale contributions show up as inverse powers of e only when one evaluates average values of 
those operators over a solution of the Schrodinger equation in d = 3 — 2e dimensions. The 



crucial observation made in Refs. [ll,12j| for the 0(ma ) corrections to positronium levels is 
that even without knowing an explicit form of this solution but using only the Schrodinger 
equation itself, one can extract all the divergent pieces in the form of (S(r)) /e, where r is 
a distance between the electron and the positron while the average value is calculated over 
the c?-dimensional wave function. Since after such an extraction the divergences contained 
in both hard and soft scale contributions have exactly the same form it is easy to check that 
they cancel each other so that a finite remainder can be safely calculated in three dimensions. 

I employ the same idea for the helium atom, where an analytic form of the wave function 
is not available even in three dimensions. Nevertheless, in perfect analogy to the positronium 
case, the use of the Schrodinger equation alone makes it possible to extract the divergent 
pieces of all soft scale contributions on the operator level. Performing such an extraction 
I manage to demonstrate straightforwardly that the divergences coming from both scales 
cancel each other before any numerical calculation. As the result, the total 0(ma 6 ) correc- 
tion to a singlet S level is represented as a sum of apparently finite average values of the 
regularization-independent operators. These average values can be calculated using a wave 
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function of the helium atom, built as a numerical solution of the Schrodinger equation |J. 

It is worth mentioning that the idea of the approach has a simple physical reason. In 
fact, soft scale divergences in bound state energy are of the ultraviolet origin. Hence they 
should be proportional to a value of the corresponding wave function at zeroth separation 
between interacting particles. In terms of the effective theory it means that by virtue of 
the Schrodinger equation one should be able to rewrite the singular soft scale contributions 
in such a way that corresponding divergences are shifted to the Wilson coefficients of the 
contact operators. After the perturbation theory is reformulated in such a manner, and if 
the underlying theory is renormalizable, all divergences that appear in any given Wilson 
coefficient have to cancel each other. 

Validity of the results obtained below for the helium can be checked in two limiting cases. 
The first, a — > at finite Za describes helium-like ion with the electron-electron interaction 
switched off. The second, Z — > at Z2 — > 1 describes parapositronium. Since in both 
cases three-dimensional wave functions of all S states are available in an analytic form, the 
average values of effective operators can be calculated explicitly (modulo (S(r)) je terms). 
Comparison with the known results shows complete agreement for all contributions. 

In order to make the formulae more transparent, I write the nonsingular soft scale oper- 
ators with coefficients taken at d = 3. 



III. ORDER mo 5 EFFECTIVE OPERATORS 

This Section illustrates the general scheme by the calculation of effective operators in 
the first non-trivial order. As previously mentioned, at 0(ma 5 ) the relevant contributions 
to the energy come from three scales: ultrasoft, soft and hard. Below we will calculate 
corresponding effective operators. 



A. Hard scale contribution 



Hard scale effects in the interaction between nonrelativistic particles a and b give rise to 
the contact operators which show up as c a bS(r a b) in the spatial representation and therefore 
as c a b m the momentum one. In the NRQED approach, c ao is extracted through the matching 
procedure, namely, equating the ab — ► ab scattering amplitude calculated in the full QED 
to that in the effective theory approach. Dimensional regularization is best suited to this 
procedure: c a b equals minus QED scattering amplitude for the particles a and b taken on 
their mass shells and at rest. In this manner the hard one-loop vertex correction to the 
single Coulomb exchange between the spin-1/2 point particle a and a particle b induces the 
effective potential^ 

V hC (r ab ) = (- -2\nm a ) 5(r ab ). (7) 

3 mi Ve / 



2 To simplify the presentation, I omit the factor (4-7r)T(l + e) from the final expressions for all 



operators in Sect. [II]] . This factor does not contribute to the finite total 0(ma ) energy correction. 



The electron's anomalous magnetic moment taken in the leading one-loop approximation 
gives rise to the potential between electrons a and h 

8a 2 

Vhm(r a b) = -tt^ s a s b S(r ab ), (8) 

where s a is the spin operator of the a th particle. If the nucleus spin is zero, the corresponding 
electron-nucleus potential vanishes. The next effective potential is generated by the hard 
one-loop box diagrams. For two spin-1/2 particles such a potential reads |13| 



t t , \ {az a z b ) 2 I 1 m a + m b - 2fi ab (l + 4s a s b ) m a \ 

V hox (r ab ) = - - ln(m a m b ) - - + In — 

m a m b \ e 3 m a -m b m b I 



S(r ab ), (9) 



where fi ab = m a m b /(m a + m b ) is the reduced mass of the pair. With the 0(m/M) precision, 
the corresponding electron-nucleus effective potential is 

Hox(r eN ) = - 2 mm - ±) 6(r cN ). (10) 

The last 0(ma 5 ) contribution coming from the hard scale appears due to vacuum polar- 
ization. In ordinary few-electron atoms, an account of the electron vacuum polarization is 
sufficient for the present-day accuracy: 

KpK.) = (11) 

The hard scale contribution to the energy equals the average value of \4c + ^4m + Vbox + K P 
summed over all pairs of particles. 



B. Ultrasoft scale contribution 

According to the standard rules of the perturbation theory, virtual transition of the 
atom into an excited state induced by the emission and subsequent absorption of a photon 
is described by the following effective operator: 

f d d q Ana^. 6^-^ 

J (2^-27 ? 3ai eXp (~^) E-H-q ^ 3 * eMtqnl (12) 

where H is the c?-dimensional Schrodinger Hamiltonian, E is its lowest eigenvalue, q is the 
photon momentum, q = \q\. Assuming that q ~ ma 2 , we can restrict our attention to the 
electric dipole transitions, i.e. replace the exponents exp(±qr a ) by 1 and the current density 
operators j a by their orbital counterparts taken in the leading nonrelativistic approximation, 
j a — > z a p a /m a , in (|T2|). Integration over directions of q then gives for the ultrasoft scale 
effective operator 

aQ d d-1 f°° dqq d - 2 
Kus " (27^-1 d Jo E-H-q J ' [L6) 
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where the operator J is defined as ^Z a z aPj m a an d &d = 2ix d l 2 /T(d/2) is the (i-dimensional 
solid angle. Analytic continuation of the integral over q from the stripe 1 < Re(d) < 2 reads 

dqq d - 2 _n(H-E) d - 2 



lo E — H — q sin(7rc?) 
Expanding now the right-hand side of ( |I3| ) in e = (3 — d)/2, we get with the O(e ) accuracy 

V ns - ^J(H -E)(± + l-\n 2(H - E)) J. (15) 

As previously mentioned, the 1/e term is due to the divergence of V ns in three dimensions. 
The 0(ma 5 ) ultrasoft scale contribution to the energy is the average value of the operator 
(|15|) over the <i-dimensional wave function: 

Here I used the Schrodinger equation and also the standard notation for the Rydberg con- 
stant Ry = ma 2 /2. Since the Poisson equation [p a , [p a ,C a b]] = 4iraz a Zb5(r a — r&) for the 
Coulomb potential C a b between two particles holds in any dimensions, we have: 

[J, [H, J]] = -4™ £ z a E ( — ~ —) 2 S (r ab ). (17) 
The ultrasoft correction (lift) in helium then reads 



where the helium Bethe logarithm WM is defined as 



k n (p 3 (H - E) In £=% 3 

In = ^ , 7 ^ , % 3/ (19) 

Ry (p 3 (H-E)p 3 ) 



and can be safely calculated in three dimensions. 



C. Soft scale contribution 

At the soft scale, momenta of all particles (electrons, nucleus and virtual photons) are 
of the order ma. Therefore the 0(ma 5 ) soft scale effective operators are generated by 
transverse photon exchange(s) since only such scattering amplitudes can contain odd powers 
of photon momentum. 
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1. Retardation 



Let us again start with the single transverse photon exchange described by the effective 
operator (|T^) but now assuming that q ~ ma. Since H — E <C ma, we can expand 
the integrand in (H — E)/q. Zeroth order term of this expansion describes the magnetic 
interaction in instantaneous approximation and is included (modulo relativistic corrections) 
into the Breit perturbation (see (|4])). The first order retardation effect is represented by the 
operator 



Vr 



ret 



Ana 



d d q &ii ~ "-f 



(2ir) d 2q* „ 
Using again the Schrodinger equation, we get 

l q S ij 



J^jai exp(-iqr a )(H - E)Y,3bj exp(iqr t 



(20) 



8 iet E = ixa 



A d n A. . _ liSi 
a a u v q 2 



(2tt) c 



Q 



^Zjai exp(-iqr a ) 



[21) 



The order ma 5 correction arises due to the nonrelativistic current densities, 



3a\P,P) 



Z a 



p' + p + 2[(p' - p)s a ,s a 
2m 



(22) 



We then see that in (21) only the exchange between different particles (a ^ b) can give a 
nonzero contribution. In fact, the integral over q in the 'diagonal' terms (a = b) is scaleless 
and hence vanishes. The expression (^TJ) simplifies, 



z a z b 



6%E = -2na £E 

a b> a m * m b 



(Uij(r ab ) \pai, \pbj,C ab (r ab )]}) 



(23) 



and reproduces the result obtained in Ref. |T3 for positronium. Here 



Uij{r) 



d d q Sij 



(2vr)< 



Q 



exp{tqr) 



d-3 \ r 3-d 



d-3 



67r (d+l)/2 



2 



-riiUj , 



(24) 



C ab (r) = az a z b T(d/2 — l)r 2 d /Tc d / 2 1 is the Coulomb potential in d dimensions and n = r/r. 
Thus we get 



4 ^r(^)r(f) 



7T 



d-3/2 



7 2 7 2 , 

Z a^b / r 3-2d 



\* <>i, 

a b> a " L a" L b 



(25) 



Here we cannot take the limit d — > 3 since the average value of r 3 diverges logarithmically. 
However, we can extract the divergence in the following way. By definition, we have 



,3-2d 
ah 



f f f 00 dp 

&at J dr' J dn ab -^^^{^pn^r'). 



/o p 



(26) 



Here a. ab = \z a z b p ab a\~ l is the Bohr radius for a given pair of particles, r ab = & ab pn ab being 
their relative position. The remaining independent variables are denoted by r'. Integrating 
by parts in the last integral, we get: 



S 



roc (j^p ]_ roo Q 

/ — r ^-ip 2 {& ab pn ab ,r') = —— dp p 2e —ijj 2 (a ab pn ab ,r'). (27) 
Jo p l ~ Zt 2e Jo op 

Here I took into account that lim^o P 2e being written as the integral over momentum has 
no scale and hence vanishes. Substituting p 2e — > 1 + 2elnp + 0(e 2 ) into the r.h.s. of (p7|) 



gives for the retardation correction (25j): 



If" 2 >» — - Wr-)> - ) , (28) 

3 „tri m a m b V e a ab J \ Ttri h 



a b>a 



where V a t = d/dr ab and 7 = 0.5772 ... is the Euler constant. The arrow shows what the 
gradient is acting on. Thus we have managed to extract the divergences in the form of 
average values of the contact operators S(r ab ) divided by e. Since if) and its first derivatives 
are finite^ for r ab — > 0, the non-contact average values in (|28|) are finite in three dimensions. 

2. Double seagull 

One more soft scale contribution of the order ma 5 appears due to the double transverse 
exchange between two particles when both photons are emitted and reabsorbed in the seagull 
vertices. The corresponding effective potential derived in [IS] for the positronium can be 
easily generalized to a more complex atom: 



V ds = -2a 2 



rn-e 



,2 



(47T) 1 " 



2c 



e 17 - 8h 'W) 



z„ z h 



a b > a '"'allt'b 



Exploiting the same trick as above to extract the divergences we get the double seagull 
contribution to the energy: 



(3) 2n „ z 2 z 2 J/1 , ll-2In2\ /P , NV /7 + lnff ^ \ 
€ ] E - -« 2 E E -f^- 7 + 2 In a ab (S(r ab )) - ( -» 



(30) 



D. Total ma correction 



The e _1 terms cancel out in the sum of all 0(ma 5 ) corrections to the energy. Hence 
we can take the limit d — > 3 in this sum. With the 0(m 2 /M 2 ) precision, the result for the 
helium ground state reads: 



3 It follows from the Schrodinger equation in d dimensions. 
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^ = -^(*Mz<.) + 2h.^-H) Wri) + , W > 

2 < Za ^ln(Z Q ) + 41„ *° 



3mM V v y Z 2 Ry 



3irmM \ r\ r\ 




7(Za) 2 1 7 + ln(Zmar 1 ) -± j + \n(Zmar 2 ) =t \ . . 

+ - — \ 3 71,1 V 1 "I ^2 n 2 V 2 / • ■ > 1 ) 



Here and below I use simplified notations: ri = ri/ii = r 2 3, T2 = r2?i2 = 7*31, r = rn = 7*12, 
the gradients are taken over the corresponding position vectors, Vi = d/dri and so on. In 
the limit of no recoil (m/M — ► 0), the result ( |3"TD agrees with the results of Araki |T(| and 



Sucher [17] after integrating by parts in their average value Q: 



« - S3 + h + HmapMr)) = -± (l±« „V) . (32) 



The first recoil (linear in m/M) correction was previously discussed in Ref. [IS 



IV. ORDER ma e HARD SCALE CONTRIBUTIONS 

Similarly to what was done in the previous order, one has to consider the hard scale part 
of a two-particle scattering amplitude but now in two loops. There is no need to consider 
t/iree-particle scattering amplitudes. In fact, the probability density to find three particles 
forming the helium atom at the same point is of the order (ma) 6 . On the other hand, these 
particles should exchange at least three photons to form a hard loop. Hence, hard scale 
effective operators proportional to S(r a b)5(rb c ) can produce an 0(ma 9 ) correction only. 

The radiative recoil potential appears when we account for the first radiative corrections 
to the hard one-loop box diagrams (see, e.g., WW). The corresponding two-loop diagrams 
involve only even powers of the electric charges Z\ and zi. Hence the radiative recoil effective 
operator coincides with that for parapositronium ||T9| . |T2f : 



^rad rec 



/6C(3) _ 1099\ 
V 7T 2 27tt 2 72 J 



na 3 
m 2 



(5(r)) . (33) 



The corresponding electron- nucleus operator vanishes in the non-recoil limit m/M — > 
considered from now on]]. Then, one- [20] and two-loop |2T|-p3[| pure radiative corrections 



to electron-nucleus interaction give rise to the following energy shifts: 

C^=(f "2 1.2)^(^0 + ^ 2 )>, (34, 



4 Order m 2 a 6 /M correction is much less than unknown ma 7 corrections. 
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re N r, I 9C(3) 2179 31n2 10\ na 2 (Za) ,_, , „ 



The net effect of two-loop contributions to the slope of the electron Dirac formfactor pT 



Pauli formfactor [22|, and vacuum polarization [23] reads: 



, ee ^ /15C(3) 631 rl n 29\ vra 3 ... .. 

Finally, to get the pure recoil contribution to the electron-electron hard scale interaction 
(three photon exchange) we have to change the sign of the corresponding parapositronium 
result^ [0: 

r ^ ( 1 . i 39 C(3) 32 „, n 7\ vra 3 ... 

S rec E= — + 41nm ^-L + — -6 In 2 + - -— 5 r . 37 

ye 7T^ 7H 3 J Am 2 

Among the hard scale contributions only the last one contains the divergence. 

V. ORDER ma 6 SOFT SCALE CONTRIBUTIONS 

The aim of this Section is to demonstrate that in analogy to the previous order the sum 
of all singular average values reduces to the form 

«^ = ' (38) 

so that the sum of soft and hard scale contributions is finite in three dimensions. 

There are many soft scale effective operators with singular average values. One can easily 
determine whether an average value of a given operator is singular or regular for d — > 3 using 
only the fact that the wave function and its first derivatives are finite when positions of two 
particles coincide. 

A. Irreducible corrections 

1. Dispersion correction 

Let us consider details of the singularities extraction procedure using the dispersion 
correction as an example. Nonrelativistic expansion of the electron's dispersion law, u p = 
y/m 2 + p 2 , reads: 

2 4 fi 

u = m + £_ - J— + p + (39) 
p 2m 8m 3 16m 5 



5 Recall that it is convenient to omit the overall factor (47r) 2f T 2 (l + e) from the final expressions 
for all 0(ma 6 ) operators. 
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The last written term induces a correction of the appropriate order. Using the Schrodinger 
equation we get 



P1+P2 
, 16m 5 



2m? 



'pI + pI\ 3 _ 3 pIpI p\ + pV 



2m 



Am? 2m 



1 1 {E - C){H - C){E - C) -3<f^±4,£-C 



2m? 
1 

2m? 



(E-cy 



c. 



[H,C] 



Am 2 



m' 



(40) 



Here the total Coulomb potential C is the sum of the electron-nucleus and electron-electron 
parts, 



C = C N + c = C x + C 2 + c, 



(41) 



while a pairwise Coulomb potential is defined after Eq . ([24]) . Singular contributions to (g0| 
are induced by the following operators: i) C 3 ; ii) the double commutator, 



2m 2 



where 



C. 



[H,C] 



2m? 



c A S a V a , eV 



a=l 



2m 



m 



{S 1 -S 2 )e £ 2 + £, 



2m 3 



+ 



Am? 



+ 



2m 3 



e = - [V, c] 



£ a = -[V a ,C a ]; 

are electric forces exerted on electrons; and iii) the anticommutator 
3 / f /?!?> ) \ 3 



(42) 
(43) 



16m 2 



' p\p\ 
m? 



' 2 1 2 \ 2 

p{ + pi\ 



' 2 2 \ 2 
p\-pi\ 



16m 2 
3 

16m 2 \ 



2m 



2m 



c 



{2c(E-C) 2 -[c, [H,C}\ 



c, 



(Ppf 



m c 



(44) 



where P = p l + p 2 and p = (p 1 —p 2 )/2. Double commutator [c, [H, C}} can be transformed 
similarly to (§3), the last term from (IP]) is conveniently rewritten as 



16m 2 



(Ppf 



m z 



16m 4 



([Pp, [Pp,c]] + 2(Pp)c(Pp)) 



(45) 



Summing up all of the above contributions and using the virial theorem in three dimensions, 
(C) = 2E, we get for the dispersion correction: 



5E 3 3E 2 (c) 



2m? 



5m* 



3E 
2m? 



, Cc v\p\ \ 



C 2 - — 
2 



Am 2 J 



, / MCnpI 3 [Pp, [Pp, § 3 (Pp)c(Pp) 9C 2 c 3C N c 2 



Smr 



16m 4 



5m 1 



Am? 



(46) 



5(£i-£ 2 )e C| 



16m 3 



2m 2 8m s 



+ 



+ g 
Am? 



+ 



5m J 
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One can easily check that only the operators C%, c 3 and £%, e 2 have divergent average values 
in three dimensions. The following analysis shows that in a similar manner singularities of 
all soft scale effective operators appear either as the third power of the Coulomb potentials 
or as the electric forces squared. 



2. Coulomb corrections 



The order p 4 /m 4 correction to the Coulomb potential between electron and nucleus reads: 



V, 



cN 
C 



128m 4 



pl\plCx] 



3nZa 
16m 4 



{pl,6(r 1 )} + (1^2). 



(47) 



I drop the spin-orbit term, which vanishes in a state with the total spin zero. The average 
value of the double commutator can be conveniently rewritten as 



pilpid] + (l<-2) 



2m([tf-C7,[p?,Ci] + (l<->2) 
-4m lE\ + E\ + [S x - £ 2 )e) . 



(48) 



For the anticommutator in (fi7), it should be taken into account that the average value of 
the operator 5{r a )C a is a scaleless integral over r a and hence vanishes in the dimensional 
regularization: 

(49) 



[{pl,5( ri )} + {pt,5(r 2 )}) = (tf( Pl ) (Am(E - C 2 - c) - 2p\) + {l^2) 
The energy shift induced by (|47|) is 



32m 3 



+ (1 - 2). 



(50) 



Note that e changes its sign under the permutation (1 «-> 2). 

Similar analysis for the correction to the electron-electron Coulomb interaction, 



C 



128m 4 



Pi,[Pi^]]+^{p?,5(r)} + (1^2) + 



(PlP 2 ) C (PlP2) -PliP2jCp lj p 2i 



32m 



16m 4 



(51) 



gives 



5%E 



' 3 [Pp,[Pp,c}} | Tmgrj 



32m 4 



— Cm — 



m- 



3P 2N 
32m 



(gi ~ g 2 )e 

16m 3 



5m J 



(52) 



Virtual transitions of electrons to negative-energy states induced by the Coulomb exchanges 
generate the energy shift 



5 C E 



'£f + Sl + 2(£ 1 -£ 2 )e + 2e 2 



5m J 



(53) 
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The total irreducible Coulomb correction is the sum of (|5UD, ([52]), and 

Ena 



5 C E 



Am 3 
1 37iaZ 



(45(r)-3Z [<*(>•,) + ,5(r 2 )]> 



1 



+3 



\ 4m 3 

[Pp,[Pp,c\] , (£i-£ 2 )e + 



Tracer) /3P 2 



m J 



32m 4 



+ 



32m 3 



32m 3 



+ 



5m 



32m 



(54) 



5. Magnetic corrections 



These corrections originate from single magnetic exchanges between particles in instan- 
taneous approximation. There are two sources of such effects. Relativistic corrections to 
the instantaneous interaction of the Pauli currents of the electrons induce the following 
contribution to the energy: 



Px+pI P\cp 2 + (d - 2)(p 1 n)c(np 2 ) iraS(r) 



2m 

d-1 
~ 32m 4 



Pi 



P 2 i,c 



4m 3 
+ 



pi, 



Pl,c 



which can be transformed to 

E I Picp 2 + {p 1 n)c{np 2 ) 2na5(r) 



5%E 



m 



+ 



/ 2%a8(i 



2m 2 



-a 



\ C p 1 cp 2 + (d- 2)(p 1 n)c(np 2 ) \ 
\ Am? ' m J 



[Pp, [Pp,c]] (e 1 -e 2 )e d-le 2 



irrr 



(55) 



(56) 



Then, virtual transitions to negative-energy states of one electron induced by a single mag- 
netic exchange with the other and a single Coulomb exchange with both other particles, 
shift the energy by 



5- M E 



£ 2 )e d-le 2 



2m 3 



2 m 3 



(57) 



4- Retardation corrections 



The effect of retardation on a single magnetic exchange between the electrons is described 
in the appropriate order by the following average value: 



SrE 



d d q 2na 
(2n) d g 4 



H e -iqn Pi + 



rn 



H, e 



qr 2 Pl ~ [ g 2g, s 2 ] 

m 



+ H.c. 



(58) 
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where q is the magnetic photon's momentum and p q a = p a — q{p a q) / 'q 2 • The correction fl58] 
consists of zero-, single-, and double-Coulomb parts: 



(2vr) rf 


m 4 


d d q 


2ira 


{2n) d 


m 3 


d d q 




(2n) d 


m 2 



5 R E = A + Ai + A 2 

, , d d q na 
A > 1 

Ax 

A 2 



' 2 PlP2 
P2i 



d-l 2 

4 g c »qr 



[CPi] 



2 P2 ;qr 

P2, g4 e 



+ (1^2) 



(59) 
(60) 

(61) 

(62) 



Integrating over q and using the Schrodinger equation we transform these expressions to the 
following ones: 



A r 



E 2 (c) , E(Cc) 
4m 2 



8m 2 



+ 



+ 



C^c Cnc 2 7ra5(r 



5m' 



4m 2 



+ 



2m 4 



-P 2 + 



pep 

4m 4 



P 2 



(Pp)c(Pp) p(c(np 2 ) 2 + (p^yep 2 , - 3(p 1 n) 2 c(np 2 ) 2 + (1^2) 



5m 4 



16m 4 



+ 



(Si-e 2 )e 



c 3 rf-le 2 
^7 + —; ; 



5m z 



m J 



(63) 



Ai 



£ (c 2 ) + / cr 2(np 2 )(£ lP2 ) + (wgi) [(np 2 ) 2 - p 2 ] + ^ 



4m 2 

(p x n) 2 c 2 + c 2 (np x ) 2 



16m 3 

:i<- 



N C N c 2 d - 2 c 3 d - 1 e 2 
2) - ; n + 



A, 



cr 



16m 3 4m 2 4 m 2 4 m 3 

,3S 1 £ 2 - (n£ 1 )(n£ 2 ) - 2{£ x - £ 2 )e (d - l)(d - 2) 2 c 3 \ 



5m z 



$(4 - d) m 2 



(64) 
(65) 



Here and below I use short-hand notations (np a ) 2 and (p a n) 2 for the operator ni(np a )p a 
and its hermitean conjugate, respectively. 



5. Seagull correction 

Double magnetic exchange between the two electrons one of which goes over into 
negative-energy intermediate states gives rise to the energy shift 

A r? / Pi c2 Pi + 3( Pl n)c 2 (n Pl ) d-l e 2 \ 

OsezgE = ( — -3 + (1 «-> 2) + — : 3 ) • 66 

This correction completes the list of the 0(ma 6 ) irreducible contributions to the helium 
energy. 
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B. Reducible Corrections 



1. Breit Hamiltonian 

Breit Hamiltonian for the helium singlet states, 

U = U S + U P , (67) 

consists of two parts, Us and Up, with the selection rules |A«S| = and |A«S| = 1, respec- 
tivelyQ. Just like in the positronium case (see |T2[]), the second order iteration of the .S-wave 
Breit perturbation diverges in three dimensions. That is why this perturbation should be 
considered in d dimensions: 

U$ = _4±|2 + ^ ( , (ri) + S (r 2 )) + Mz3l5, (r) - P^ + { d- 2)( Pin) c( np ,) 

8m d 2m z m l 2m z 

(68) 

As for the P-wave part, which mixes singlet S and triplet P states, the corresponding second 
order iteration is saturated by the soft scale and therefore we can take this perturbation in 
the limit of d = 3: 

«i - « 2 (Zh Zl 2 rx P \ 

Here l a = r a x p a . Below I consider energy shifts arising in second order in Up and Us- 

2. Second iteration of the Breit Hamiltonian: P-wave 

To find the singlet S level shift induced by an admixture of triplet P states, 

5 P E = (UpGUp) , (70) 

where G is the (reduced) Green function of the Schrodinger equation, consider first the 
action of Up on the ground state wave function. As far as the latter depends on the absolute 
values of r 1; r 2 and r, tp — ip(r 1 ,r 2 ,r), we can substitute: 

li = r± x Px — > —ir\ x \ri\0\ + no) = i o, l 2 — > —i o, 



r x P — > — zri x r 2 ( — <9i H <9 2 

ri r 2 



t/p - ia A ^- 82) (4-d - + (1 - 2)) , (71) 

2m z yrfr r^ri I 



where d a = d/dr a , d = d/dr, while the vector 



si + s 2 



16 



ri x r 2 r x ri r x r 2 



(72) 



is perpendicular to the triangle composed by r±, r 2 and r, while its norm equals this triangle's 
area. The P wave admixture to the S state wave function, 



5 p ip = GUpip, 

has the same angular dependence as the perturbation: 

A(s 1 - s 2 ) (j)p(r 1 ,r 2 ,r) 



5ptp = iot- 



2m? 



r\r 2 r 



(73) 



(74) 



Function </>p(ri, r 2 , r) is introduced in such a way in order to make it finite at coalescense 
points: 4>p(ri,r 2 ,r) < oo when r a — > or r — > 0. Substituting the r.h.s. of Eq.([74|) into the 
inhomogeneous equation, 

(E - H)5 P 4) = Upij, (75) 
we can cancel the spin- and angular-dependent factor using the following relations: 
[H, 8l - a 2 ] = 0; 



[H, r 1 x r 2 ] 



_Pi_ 
2m 



ri 



x r 2 + ri x 



2m 



r 2 



-<9i H — <9 2 + -a ) . 



m 



At the last step I took into account that the operator acts on the function depending only 
on r 1; r 2 and r. The resulting equation for <pp reads: 



\H> --(d 1 - + d 2 - + d-) -e\<Pp = \ T -\d, 



2 T d 

r 2 



2)U. 



(76) 



Here H' is obtained from H by the substitutions d — ► d—l/r, di — > d\— 1/rx, d 2 — > <9 2 — l/r 2 . 
The final expression for the energy shift (|70|) then reads: 



5pE 



a 



4m 4 



A 2 



^ 2 -^+(l-2) 



(77) 



5. Second iteration of the Breit Hamiltonian: S-wave 

In order to extract divergences from the second order correction induced by the .S-wave 
part of the Breit Hamiltonian, Eq. fl6"8|) , it is convenient to rewrite the latter in the form 

U s = + {H, u} - 2Eu + O d , (78) 

where 
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M=q , + («l-l) e 

4m 

q, + (J-i )e + _i_f ! i + jg + 1_ gg£_2^3 c2 

2m 8m [ m m J 2m 4m 
p x (C N - (d- 2)c)p 1 + (1 ^ 2) p l cp 2 + (d-2)(p 1 n)c(np 2 ) 
+ 4m^ 2m^ ' (8 ° } 

Inserting ( [78]) into equation 

5 5 £ = (C/5GC/5} , (81) 

and using the Schrodinger equation for the reduced Green function, (H — E)G(R, R') = 
Tp(R)ip(R') — 5(R — R'), where R denotes the vector (ri,r 2 ,r), we get: 

S S E = 2 (U S ) (u) - (~ [[H,u]M + {U 8 ,u}) + (OsG0 3 ) . (82) 

Note that we can take the limit d — > 3 in the last term since in contrast to Us the perturbation 
O3 does not contain operators more singular than C 2 and c 2 . In other words, all the 
divergences in (|8~2"D are moved to the average values of local operators. 
The first term from Eq.(|8"2"D is most easy to calculate: 

2 (U 8 ) (O,) = B 2 -^±^. (83) 

Here B = (Us) is the non-recoil limit of the first order Breit correction (|j). For the second 
term from Eq . (|82|) we have 

/Irrrr 1 A /£i 2 + £f (S l -£ 2 )e (d-l) 2 e 2 \ . x 

Then, the third term from Eq.(^) can be rewritten as 

/ 3CcC N p 2 C N p 2 [Pp,[Pp,c}}+2(Pp)c(Pp) 
\ Am 2 8m 4 8m 4 
ira6(r) c | j C N + (d-l)c p 1 cp 2 + (d - 2)(p 1 n)c(np 2 ) \ 
2m 3 [ 8m 2 m J 

(S 1 -£ 2 )e C 3 d-1 3 E 2 + g (rf-l)e 2 \ 



4m 3 4m 2 8m 2 8m 3 8m 3 



Finally, the fourth term from Eq. (|82|) can be calculated in a way similar to that used above 
for the P-wave contribution. Namely, we first find the solution for the inhomogeneous 
equation 

{E-H)8 s ip = {0 3 -(0 3 ))iP, 
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orthogonal to the ground state, (5sip\ip) = 0, and then evaluate the matrix element of 3 
between 8 sip and ip: 

(0 3 G0 3 ) = (8 s ip\0 3 \iP). (86) 

An alternative way to calculate (]86f) is to include the operator 3 directly into the Hamil- 
tonian: 

{H + 3 W = (87) 

Then, with the 0(a 2 ) precision, 

(0 3 G0 3 ) = E' — E — (0 3 ) . (88) 



C. Total soft scale contribution 



Summing up all soft scale contributions, given by Eqs.fl4q, |54| , pq , p7| , |63H66|, \n\ , 

we get: 



8 so ftE 



E ? - 



+ 



2m. 2 
B{c) 



E 2 (c) , £ /2 C N C + c 2 ttqZ 
4m 8m 3 Am 2 



Am 2 ^ m 



[8(ri) + 6(r 2 )} 



2m 



5 P E + (0 3 G0 3 ) 



cC N C f C N p^Pa + (p 1 n)c(np 2 ) 1 p?C N p| 



2m? 



4m 4 



6m° 



Am A 



p?c(np 2 ) 2 + {p^nycpl - 3(p 1 n) 2 c(np 2 



or 



16m 3 



16m 4 

2{np 2 ){S 1 p 2 ) + (n£ 1 ) \{np 2 ) 2 - p 2 2 ] +(1^2) + H.c. 



cr 



(3£ 1 £ 2 - {n£ 1 ){n£ 2 ) - 2{S 1 - £ 2 )e) 



8m' 

3 [Pp, [Pp, c}} + na5(r) ( 13P 2 + ^ 



32m 4 



2m 3 



16m 



Am 3 

r 3 



4m 2 



+ 



5(ri 



3p& 



— ^ + 2C 2 + c + (1 <-> 2) 
2m / 



(gi-g 2 )e 

32m 3 



1 — 3e 3 3 — 6e 2 



5m J 



2m 2 



Am 3 



(89) 



Order e coefficients are kept in Eq.(^) only if they multiply the operators whose average 
values contain 1/e singularities. In (|89|), the following relations are taken into account: 
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[P., fe, Ml = ( d - 3) ({c /^^y^ }) = (i - 3) (c=) + 0(e). 

(90) 

For the bulk of the operators in fl55|), their average values can be safely evaluated at — > 3. 
Special care is needed when one deals with the operator p 2 c(np 2 ) 2 + (1 «-> 2) + H.c. which 
is not well defined in three dimensions. We can take the limit — ^ 3 having previously 
determined what the momenta operators are acting on: 



lim (plc(np 2 ) 2 + (1 2) + H.c.) = (pjc(np 2 ) 2 + (1 2) + H.c. + ATrm 2 a 3 5(r 

d ~>3 \ ' \ ' I d=3 

(91) 

In order to be certain that this relation holds, it is sufficient to consider the two-body 
problem, where both sides of Eq. can be calculated analytically. 

In order to calculate singular average values entering into ([59]), consider the Coulomb 
potential C ab between two particles with charges z a and z b . The average value of C 3 b is 

'A) = (W0 3 ( T -^S\ (r 6 ab 3d )- (92) 



7T 



Repeating the procedure used in Sec. |TJ we get 

(C 3 ab ) = {z a z b af | Q - 41n^- + 2) (nS(r ab )) - 2 ( j + ^ ^ n a tf ab ^\ . (93) 

Recall that a. a b = ^aZb^ab^ 1 ■ Average value of the second singular operator, £ 2 b , can be 
found in the following way: 

(£ab) = ([V , C ab f) = ([V„, C ab [V a , Cab]] - C ab [V a , [V , C ab ]]) = -2 (C ab [V a , C ab ]^ a ) . 

(94) 

At the last step, I again have used the equation (C a bS(r ab )) = valid in the dimensional 
regularization. In order to express the average value 



£ 2 ab ) = 2(d-2)(z a z b a) 2 [ T [l 2 _ 1 l J ) (r*- 2d n a rf ab ) (95) 



7T 



in terms of (Cf 6 ), the known derivative of the wave function at r ab — > 0[] can be added to 
and subtracted from n ab ^ ab : 

S 2 b ) = 2fi ab (d - 2) (c^) + 2(z a z b a) 2 ^^J- (n ab ^ ab - n ab z a z b aj ^ . (96) 



Since the last average value in the r.h.s. of fl96|) is finite, it can be considered in three 
dimensions. Now, extracting all divergent pieces from ([59]) we get (RHJ). 



7 Strictly speaking, with a coefficient 1 + 0(e). Inspection of the known two-body average values, 
however, shows that the coefficient is in fact equal to 1. 
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VI. CONCLUSION 



The divergences contained in the hard (^?p and soft 
other so that in the sum of all contributions, Eqs. (IST 



account Eqs. (|9ll , |93| , |96D , the final expression for the 0{ma 
energy of the helium atom can be written as 



scale contributions cancel each 
i, we can put d = 3. Taking into 
) correction to a singlet S-state 



E 3 



2m? 
B(c) 



+ 



E 2 (c) 
Am? 




2C N C + c 2 
4m 



PiPl 
8m 3 



traZ 
4m 2 



2m 



+ 5 P E + (0 3 G0 3 ) + vraW (k eN (5( ri ) + 5(r 2 )) 




+ 



3CiC 2 Cn cC n C C N c[p 1 p 2 + n(npi)p 2 ] + H.c. 



4m 2 



2m 2 



irrr 



pICkpI , Pic 2 p 1 +p 2 c 2 p 2 , (Pi x p 2 )c(Pi x P 2 



4m 4 



8m 3 



+ 



plc(np 2 ) 2 + (p 1 n) 2 cp§ 



4m 4 

3(p 1 n) 2 c(np 2 ) 2 - 



-a- 



16m 4 

2(np 2 )(£ip 2 ) + (n£i) [(np 2 ) 2 



p 2 ] + (!<-> 2) + H.c. 



3£i£ 2 — (n£i)(n£ 2 



16m 3 

-2(£!-£ 2 )e 



■5m 



7ra5(r) /9P 2 



2m 3 



16m 



£ 2 )e + (Za) 2 



-a: 



32m 3 4m 3 
3 ln(mar) + 7 



2m 2 r 2 



nV + 




+ (1 - 2) 



2m 3 r 3 



(97) 



Here all momentum operators standing to the right (left) of position-dependent operators 
are assumed to act on the right (left) wave function. Although d- dimensional notations for 
the Coulomb potentials and electric forces are kept in (|97D, the immediate three-dimensional 
counterparts are implied for all operators, e.g., C\ — > —Za/ri, e — > an/r 2 and so on. The 
contact terms enter into Eq. (p7|) with the coefficients 

10 Z 9ZC(3) 



&eN 
h 



Z 3 A27Z 2 
~2 + 96 



In a + 



3385 



27 
331 



2179Z 3Z - AZ 2 , 

n H m 2, 

4tt 2 648tt 2 2 

29 In 2 15C(3) 



+ 



(98) 
(99) 



216 54vr 2 2 4tt 2 
Eqs. (p7j - p9|) are the principal result of the present work. Its application to the ground state 
of the helium atom is considered elsewhere 0. 

The approach elaborated here can be applied to other few-electron atoms as well as 
to higher order corrections. The only stumbling block to higher order calculations is the 
yet unknown three-loop hard scale electron-electron potential. Nevertheless, the order ma 7 
corrections enhanced by powers of In a can be determined by combination of methods used 
in Ref. fl5j and in the present work. 
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